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ABSTRACT 

We compute the covariance expected between the spherical harmonic coeffi- 
cients aim of the cosmic microwave temperature anisotropy if the universe had 
a compact topology. For fundamental cell size smaller than the distance to the 
decoupling surface, off-diagonal components carry more information than the di- 
agonal components (the power spectrum). We use a maximum likelihood analysis 
to compare the Wilkinson Microwave Anisotropy Probe first-year data to models 
with a cubic topology. The data are compatible with finite flat topologies with 
fundamental domain L > 1.2 times the distance to the decoupling surface at 
95% confidence. The WMAP data show reduced power at the quadrupole and 
octopole, but do not show the correlations expected for a compact topology and 
are indistinguishable from infinite models. 

Subject headings: cosmic microwave background, cosmology: observations 



1. INTRODUCTION 

The simplest model for the universe is a spatially homogeneous, isotropic spacetime with 
a Euclidian (flat) geometry. This simple model is consistent with observations, but leaves 
unaddressed the question of topology or the connectedness of spacetime. Schwartzschild 
(Schwartzschild 1900) first noted the possiblitity of non-trivial topology for the Universe 
even before Einstein's discovery of his field equations. Almost immediately after Einstein's 
discovery, de Sitter (deSitter 1917) pointed out the the field equations did not constrain the 
topology. Since general relativity provides no theoretical guidance, we turn to observations 
to constrain topology. 
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Observational tests of topology all rely on multiple imaging of distant objects. If the 
universe is multiply-connected with cell size smaller than the distance to some object, photons 
from that object can reach the observer via multiple paths. Simply searching the sky for 
multiply-imaged point sources, e.g. quasars, is problematic: since the travel time to each 
image is different, each image shows the same object at a different time. If source evolution is 
important, the multiple images may no longer be recognizable as such. The ideal source for 
topological tests would fill the whole sky with a pattern centered on the observer and emitted 
at a single time. The cosmic microwave background (CMB) is an excellent approximation 
to this ideal source. 

A number of authors have used the CMB to constrain the topology of the universe. These 
tests fall into two general categories. A compact topology can not support spatial structure 
with wavelength longer than the cell size. The CMB temperature anisotropy will thus be 
suppressed on angular scales larger than the (projected) cell size. The first category tests 
use the CMB power spectrum (or its Legendre transform, the 2-point correlation function) 
to test for non-trivial topology. 

The CMB in fact shows significantly less power in the quadrupole and octopole than 
would be expected for a model based on higher-order moments. The discrepancy was first 
detected by the Cosmic Background Explorer (Bennett et. al. 1996) and verified at much 
high signal to noise ratio by the Wilkinson Microwave Anisotropy Probe (WMAP) (Bennett 
et. al. 2003). Figure 1 shows the angular power spectrum of the WMAP first-year data 
compared to the best-fit ACDM model (Spcrgcl et. al. 2003). Models with compact 
topology L ~ 1 provide a good match to the observed power spectrum, motivating tests of 
finite-universe models (de Ofiveira-Costa & Smoot 1995; Luminet et. al. 2003). 

The suppression of power on large angular scales is a necessary but not sufficient condi- 
tion for the existence of a compact topology. The power spectrum is rotationally invariant, 
averaging over any phase information in the pattern of CMB anisotropy. Such phase infor- 
mation must exist for compact topologies, and forms the basis for a second class of tests. 
A "circles on the sky" search (Cornish et. al. 1998) provides a more stringent test for 
compact topologies. The CMB decoupling surface is a sphere centered on the observer. If 
the cell size is smaller than the distance to the decoupling surface, the multiple images of 
this sphere induced by a compact topology will intersect to produce patterns that match 
along circles. Such circles are not observed, limiting the cell size L > 1.7 for a wide class of 
models (Cornish et. al. 2003). 

Additional tests are possible. Compact topologies will not produce circles on the sky 
if the cell size is larger than the distance to the source, since the resulting images will not 
intersect. Compact topologies with L > 2 may still be distinguished using phase information. 
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In this paper, we describe the correlations imposed on the microwave background by the 
topology. We use this formalism to compare the WMAP first-year data to a model with 
cubic topology and derive constraints on the cell size L. 



2. COVARIANCE OF SPHERICAL HARMONIC COEFFICIENTS 

On large scales, the CMB temperature anisotropy is given by 

AT(x) = -i$(Arx), (1) 

where $(x) is the gravitational potential and At is the radius of the decouphng surface. 
The potential has the harmonic expansion 

$(x) = J ci/i(k) $k e-''^-^ (2) 

where d//(k) refiects the density of states, determined by the topology. For an infinite 
universe, k is continuous, while for a compact topology only discrete values 

k=— n, n={n^,ny,n^) (3) 

are allowed, where na;,nj^, and are integers and L is the cell size in units of the conformal 
time to the decoupling surface (Zel'dovich 1973; Fang & Mo 1987; Sokolov 1993). The 
cutoff in the discrete spectrum, |k| > suppresses power on large angular scales. 

In inflationary cosmologies, the gravitational potential's $k are random Gaussian vari- 
ables with zero mean and the covariance 

{^^K')-^S'ik-k')V{k). (4) 
We expand the corresponding temperature fluctuations 

Ar(x) = Yl «^-^^m(x) (5) 

where 

aim = J diiik) $k JeikAr) F,„(k) (6) 

and je{kAT) is a Bessel function of order £. The a^^'s have zero mean. Using Eqn (4), we 
find their covariance 

^em,e'm' = {aemal/jn') C^) 

= J d^{k) je{kAr)MkAT) ^ r,^(k)r;^,(k), (8) 
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which is not necessarily diagonal. 

In the limit of a flat open Universe, dn(k) — > k'^dkdk and the orthonormality of the 
spherical harmonics yields 

3271""^ r°° V(k) 

lim M[^^,,^, = -—5u'5mm' / dkjeikArf (9) 

For a compact topology, the integral over the continuous variables k becomes a sum over 
discrete k„ = (27r/L)n with n = {nx,ny,nz) a triplet of integers. Just as we break the 
continuous integration c?/i(k) into magnitude and angular parts, we do the same for the 
discrete case: (i/x(k) = J2neAf Xl{|n|=n} where J\f is the set of all possible magnitudes of the 
integer triplets n and {|n| = n} are all n with magnitude n. Thus '^^eAf ^^^^ 
magnitude and is, for each magnitude, the angular sum. For example, the first 

value in jV is n = 1 with the corresponding {|n| = 1} containing the six vectors (±1, 0, 0), 
(0, ±1, 0) and (0, 0, ±1). The next value in J\f is n — \pl and {|n| = v^} contains 12 vectors, 
each with ihl in two places and in the third. Writing the measure d//(k) in this form, for 
compact topologies we have 

^L,e'm' = E MknMMknAr) ^L,£'m'(^) (10) 

where 

At^,^,{n) = {-zY/ yem{n)Y,1^,{h) (11) 

{|n|=ji} 

The momentum kn is related to n via kn = xl^l- 

There is a great computational advantage for making this split in the sum. All the 
cosmology is contained in the magnitude sum while the topology is reflected in the sum 
for Afm,e'm'{n), which requires the lion's share of CPU time. When numerically evaluating 
Eqn (13), we compute and store all the required matrices Aijn,e'm'{n), which only depend on 
the relative ratios of the fundamental domain sizes L^, Ly and L^. For this work, we have 
assumed they are all the same and equal to L. As we vary the rest of the model parameters, 
including the topology scale L, we need only evaluate the factors in the first sum and use 
the stored matrices A(.rn,e'm'i'n) to complete the sum. 

Thus far, we have ignored contributions from the time evolution of the gravitational 
potential (the integrated Sachs- Wolfe effect) on large scales and from the Boltzmann physics 
of the coupled photon-baryon fiuid. on small scales. Both effects may readily be included. 
The "line-of-sight" approach (Seljak & Zaldarriaga 1996) to computing the anisotropics for 
a flat open Universe flrst calculates the (scalar) response function A^^(/i;, At) and then gives 
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the (diagonal) covariance 

MZZ''^ - i^^r^ee'Smm' fdk Ar)^ ^ (12) 

This suggests we can include the ISW and acoustic effects of the plasma motion by writing 
covariance (8) as 

ML,i'm' = (47r)^ J2 ^rhkn, Ar) Ag(A;„, Ar) ^ A^^^,^,{n) (13) 

neA/" 

with A^„,.^,„,,("') given by Eqn (11). We obtain Af^^(A;, Ar) from CMBFAST (Seljak & Zal- 
darriaga 1996), by way of the program CMBEASY (Doran 2003), along with the associated 
value of Ar. 

Equation 13 gives the covariance of the temperature coefficients function of 

cosmological parameters and global topology. Figure 2 compares the diagonal elements to 
the ofT-diagonal elements for the case of a cubic topology. The total power on the diagonal, 
Yltm \^'tmim\i reflects the content of power spectrum C/' and is rotationally invariant. The 
off-diagonal power, |, measures the correlations between the different 

angular scales and is generated by the global properties of the topology. For cell size L 
less than twice the distance to the decoupling surface, the off-diagonal elements dominate. 
Analyses based solely on the power spectrum thus ignore the main information content of 
the map. The off-diagonal correlations persist at larger cell size, but become increasingly less 
important. In contrast to the circles on the sky test, the off-diagonal correlations smoothly 
decrease past L > 2. 



2.1. Sky Map Generation 

The covariance matrix M^^^,^, for the aem coefficients fully describes a cosmological 
model and topology. Since the a^m's are still Gaussian variables, albeit now correlated, 
all higher moments are given in terms of their covariance. Equation 11 thus allows the 
rapid generation of simulated sky maps for a given topology. For a given topology scale 
L and cosmology Arp^ {k), we first compute the Cholesky decomposition ^tmi'm' of the aim 
covariance matrix: 'M.^mi'm' — ^eme'm" ^em'e'm"- Then if are a set of uncorrelated 
unit variance Gaussian variables, we set aem — X^^'m' ^imem' ^em' to obtain the sky map 
AT'(x) = J2em'^i'rn^im{'^)- A "circles in the sky" test (Gornish et. al. 1998) run on the 
resulting sky maps verifies that the algorithm correctly reproduces all features for compact 
topologies. 



-6- 



The only remaining free parameter in our model is the overall amplitude of the fluc- 
tuations. We fix this by comparing the power spectrum for the compact model to that for 
an infinite Universe with the same response function Aip^{k). The CMB power spectrum is 
given by 

= ^ E = ^ E ^Lem- (14) 

m m 

The covariance of the power spectrum becomes 

2 t I' 

^ (2£+l)(2f +1) E E ^ l^imi'm'l , (15) 

Wc normalize the covariance matrices M^^„^,„j, using the amplitude that minimizes X]^™^ {^e ~ Cf^^^^^^ 
where we use iiower = 20, as this is a small enough scale to be beyond where we expect to 
find interesting topological effects. The simulated maps correctly show the rise to the first 
Doppler peak. Since topological effects are most apparent on large scales, we take imax = 30 
for all analyses in this paper. 



3. LIKELIHOOD ANALYSIS 

We compare data from the WMAP first-year sky maps to a set of models with a cubic 
fundamental domains described by the cell size L in units of the conformal time to the 
decoupling surface. We specialize to the case of a cubic topology, identifying opposing faces 
of the unit cell without twists or rotation, so that the topology is fully specified by the cell 
size L. For each value L we compute the hkelihood log>C — —\ (x^ + logdetM) where 

= Huimm' 0-tm (M ^)e.mem' Cli'm' (16) 

and Mimffm' is given by Equation 13. 

We use temperature coefficients a£m derived from the internal linear combination (ILC) 
map from the WMAP first-year data release (Bennett et. al. 2003b). This map reduces 
foreground emission at the cost of a complicated window function and instrument noise on 
angular scales 9 < 2°. Topology is important only at much larger angular scales. We hmit 
the likehhood calculation to 2 < ^ < 30. On these scales, the dominant uncertainty is cosmic 
variance; effects from the instrument noise and beam profiles are negligible. 

The ILC map reduces foreground emission but does not eliminate it completely. We 
impose a cut in Galactic latitude |6| > 5° and compute the a£m using unit weight for all 
pixels outside the cut. Our results are stable as the cut is varied from 2?5 to 15°. 
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The likelihood analysis includes off-diagonal correlations between different aim, and is 
not rotationally invariant. Equation 13 was derived for the case when the faces of the 
fundamental domain align with the data coordinate system. Wc must thus consider different 
possible orientations between the data and the unit cell of the topology. It does not matter 
whether we rotate the model or the data, so we rotate the data. In terms of the Euler angles 
^ = (0, 9, ip), we use the rotation matrix 



where (i^j„j; (6) arc the Wigner rotation functions. Wc apply the rotation on a grid of TBD 
possible orientations ^ over the range < (}),9,ip < 7i/2. (p and 9 are uniformly distributed on 
the sphere while the azimuthal rotation ip is uniformly distributed over its range. Since the 
fundamental domain for the topology is cubic, the noted range covers all possible orientations. 

For a fixed cosmological model, the likelihood becomes a function of four parameters: 



one parameter specifying the size of the unit cell and three parameters specifying the ori- 
entation, where for clarity we have suppressed the im subscripts. Even in the absence of 
topological effects, chance alignments between the random aim can create non-zero correla- 
tions in a single realization. If the relative orientation of the data and model is held constant, 
the likelihood will approximate a normal distribution over an ensemble of CMB maps. Select- 
ing the maximum likelihood over many possible orientations of a single CMB map, however, 
selects the lowest each cell size L. If the data do not represent a model with compact 

— * 

topology, maximizing the likelihood over the (now nuisance) rotational parameters ^ will 
select from the tail of the distribution (Eq. 16), leading to a biased estimate for L. The 
bias is unimportant for large L where the scale is too large for the topology to break global 
isotropy, but becomes significant for L ~ 2 where the off-diagonal terms begin to dominate 



We quantify the effect of maximizing the likelihood over orientation using Monte Carlo 
simulations. We generate 1000 realizations of compact topologies at 51 different values of 
L on a grid from L = 0.8Ar to 4.0Ar uniformly spaced in 1/L. For simplicity, we work 
directly with the a^'m' coefficients to avoid the intermediate steps of generating sky maps, 
masking pixels near the Galactic plane, and computing the ae>m' from the un-masked pixels. 
We instead impose the Galactic cut using the matrix 




(17) 




(18) 



(Fig. 2). 




(19) 
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which removes all power from the azimuthally symmetric region \9 — t:/2\ < b, where Pp{y) 
are the associated Legendre functions. Each realization thus generates a set of correlated 

= PRL^x (20) 

where is the Cholesky decomposition of the aim covariance matrix (Eq. 13) and x is a 
vector of zero mean, unit variance Gaussian random variables. We select a grid uniform in 
1/L because the topology effects the covariance matrix M as A;„ = ^|n| oc 1/L. There is 
also the advantage of placing the infinite case at a finite distance from the region of interest, 
L ~ 1. A model with L = 3 At is nearly indistinguishable from L = oo. We vary the 
topology scale L while keeping the "background" cosmology fixed. Since we only consider 
large angular scales, the results are insensitive to the cosmological parameters. We run the 
kn integration in Eqn. (13) for n < 90, which is sufficient for convergence. 

Figure 3 shows the likelihood of the WMAP first-year data as a function of domain size 
L. When maximized over orientation, the hkelihood peaks sharply at L ~ 2. We test for the 
significance of this peak by comparing the WMAP results to 1000 Monte Carlo simulations 
of an infinite fiat universe. The simulations also peak at L = 2, demonstrating the bias 
incurred when maximizing over orientation. The WMAP data fall near the mean of the 
simulations, suggesting that the data are consistent with an infinite universe. In this case, 
the orientation ^ becomes a nuisance parameter. We may then marginalize over the nuisance 
parameter by averaging the likelihood over all orientations (as opposed to selecting the best 
orientation). The marginalized likelihood shows a plateau for L > 2.1 with a sharp drop at 
smaller cell size. Note that for L > 3 the orientation ceases to be important as the maximum 
likelihood asymptotically approaches the marginalized likelihood. 

Figure 3 shows that the WMAP likelihood falls near the mean of Monte Carlo simula- 
tions drawn from a parent population with an infinite topology. We use additional Monte 
Carlo simulations to set upper limits to the allowed size L of the fundamental domain. We 
repeat the maximum hkelihood analysis for 19000 Monte Carlo simulations, 1000 at each 
of 19 values for L ranging from 1.01 to 3.7 uniformly spaced in 1/L. For each simulation, 
we get a "best fit" topology scale Lout, the value that maximizes >C(a^'"|L, ^) (Eq. 18). We 
then tabulate the probability V{Lout\Lin) for each input realization with the "true" topology 
scale Lin to produce the best-fit output scale Lout- We invert this relationship to derive 
the probability for a given best-fit output Lout to be drawn from a parent population with 
topology scale Lj„: 

'PiLinlLout) — roc ^/ J- ry \~T77 \^ (21) 

Jo ''^{■^outl-i^inJ'J'K-i^in) 

where the factor d{Lin) explicitly accounts for the simulations' uniform distribution in 1/L. 
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Figure 4 shows both distributions. For Lin < 2, Lout — L^n and our hkehhood analysis 
successfully identifies the topology scale. On these scales, the off diagonal components 
of Mx, become dominant (Fig. 2); the correlations between the different agm values are 
important and the likelihood function strongly discriminates models. A minor identification 
Lout — Lin is also apparent. This aliasing of scale is typical and in this case corresponds 
to the ratio in sizes between a circle that just circumscribes a square and the circle that 
is just contained in a square. Such aliasing of topology scales is expected and is also seen 
during a circles-in-the-sky analysis. 

For Lin > 2, maximizing over orientation produces output Lout ~ 2.2 independent of 
the actual value of Lin > 2. The correlations between the agm induced by topology are weak 
on scales L > 2 compared to chance alignments. Individual realizations with L > 2 will 
thus appear to have a slight preference for orientation, although the parent population is 
nearly indistinguishable from the infinite model. We account for this bias using Eq. 21. An 
observed best-fit value Lout ~ 2 does not mean the parent population must necessarily have 
Lin — 2, but rather that the parent population has nearly uniform probability to represent 
any scale Lj„ > 2. 

The maximum likelihood for the WMAP data occurs for L^^^^ — 2.1. From Eqn. 

(21), the cumulative probability is P(Lj„ < Lwmapm) = Jo'" 'P{Lin\LwMAP,out)d{L-J. We 
obtain the 68% confidence limit that the topology scale is greater than 2.1 times the distance 
to the decoupling surface and 95% confidence it is greater than 1.2 times the distance. We 
place a 68% (95%) confidence that the topology scale is greater than 29 Gpc (17 Gpc). 

Figure 4a shows that a certain number of realizations with L < 2 "scatter" into the region 
Lout ~ 2 where the WMAP data show maximum likelihood. We test the null hypothesis 
(that the universe is infinite) by taking all simulations with Lout > 1-96 (the horizontal 
band in Fig 4a) and computing the cumulative probability for these simulations only. Figure 
5 shows the resulting curve. The cumulative probability for the WMAP data, computed 
using all simulations, has identical confidence intervals as the probability derived using only 
those simulations whose hkelihood peak occurs at Lgut > 1-96. We therefore accept the null 
hypothesis to conclude that the WMAP data are consistent with an infinite universe. 



4. Discussion 

A compact topology imposes a specific pattern of correlations (a^mO^'m') between the 
spherical harmonic expansion of the CMB temperature. We compute the expected corre- 
lations for the simplest non-trivial topology, the cubic torus, and compare a range of cell 
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size L to the WMAP first-year data using a maximum-likelihood algorithm. The covariance 
matrix explicitly includes the contribution from the integrated Sachs- Wolfe effect on large 
angular scales and the acoustic peaks at small scales. We separate the covariance into a piece 
dependent on the topology and a piece dependent on the cosmology. Although the transfer 
functions A^/(A;) for the cosmology assume isotropy in A;-space, which is no longer exact for 
compact models, the effect is predominantly in the cosmology with negligible effect on the 
topology. 

The algorithm is sensitive both to the power spectrum of the data (diagonal elements 
of the covariance matrix for different aim) as well as the phase information contained in 
the off-diagonal elements. For cell size L < 2 the off-diagonal elements are larger than the 
diagonal elements. A comparison of the data to topological models that utilizes only the 
power spectrum can produce false positives by ignoring the additional information in the off- 
diagonal elements. We demonstrate this using Monte Carlo simulations. The power spectrum 
is rotationally invariant and does not specify orientation. We may thus modify Eq. 16 to use 
the power spectrum Ce and its covariance (Eqs. 14 and 15) in place of the spherical harmonic 
coefficients aem- When only considering the power spectrum, the maximum likelihood for 
the WMAP data occurs at L = l.lAr; this is the "finite" model power spectrum plotted in 
Fig. 1. Does this imply a positive detection of finite topology? To test this, we generate 
1000 Monte Carlo realizations drawn from a parent population with L — 1.1 and generate 
the likehhood for each reahzation using the full covariance matrix (Eqs. 13 and 19). For such 
a small topology scale, almost all reahzations have their hkehhood peak at L — 1.1. This 
scale is small enough that the bias from maximizing over orientation is not important. For 
each realization, we also generate a "companion" realization with exactly the same power 
spectrum, but with completely imcorrelated a^m's. The two realizations by construction 
must give the same results for an analysis based solely on the power spectrum. When we 
analyze the "companion" realizations using the full covariance matrix, we obtain results 
similar to the infinite models displayed in Fig. 3. A likelihood analysis using the full 
covariance matrix successfully distinguishes models with compact topology from models with 
identical power spectra but without the correlations between different a^m required by the 
topology. Suppression of power in the quadrupole and octopole moments is a necessary but 
not sufficient condition for a compact topology. The WMAP data show reduced power at 
i = 2 and 3, but do not show the correlations expected for a compact topology and are 
indistinguishable from infinite models. 

For cell size comparable to the distance to the decoupling surface, the correlations 
become weaker. Maximizing the likelihood over orientation then allows chance alignments to 
introduce a bias in the likelihood estimator. We quantify this using Monte Carlo simulations. 
The WMAP first-year data are consistent with input models drawn from parent populations 
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with infinite fundamental domain. We establish 95% confidence hmit L > 17 Gps for the 
cell size of a cubic topology, in agreement with the result of 24 Gpc obtained by Cornish et. 
al. (2003). 

We thank G. Hinshaw for helpful discussions. NGP thanks M. Doran for his help in using 
CMBEASY. This work was supported by the National Aeronautics and Space Administration 
under the Astrophysics Data program of the Office of Space Science. 
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Fig. 1. — Power spectrum of the WMAP first-year data, compared to a model for a flat 
universe with an finite topology. The topology has a cubic fundamental domain with a side 
length L = 1.1 x the distance to the decoupling surface, the best-fit value if only the power 
spectrum is considered. The dash-dotted line shows the best-fit power spectrum for an open 
fiat Universe. The lowest three multipoles are suppressed for the finite topology model, since 
the universe in such a model is too small to support power at such large scales. 
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Fig. 2. — Comparison of diagonal vs off-diagonal power in the spherical expansion coefficient 
correlation matrix |M^^^^,^, | = | (a£r„a|,^,) | as a function topology scale L. The power 
measures the relative importance of terms. The total power in the diagonal components, 
|]^^m£ml' measures the information considered by the power spectrum and is flat with 
only a small decrease at small L due to suppression of power at the largest scales for small 
topologies. The off-diagonal power, Yli{im)i^{t'm') l-'^^m^'m'l' ^I'ises from the global structure 
present for finite topologies. This varies from almost negligible at the largest topology scales 
to being the dominant contribution for small topologies. An analysis looking for topology via 
a spherical expansion must consider the correlations between different expansion coefficients. 
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Topology Scale 



Fig. 3. — Likelihood of the WMAP first-year data as a function of topology scale. At each 
topology scale L, the likelihood is maximized over the possible orientations of the funda- 
mental domain. For L ~ 2, chance alignments create a bias in the likelihood estimator. 
The solid line and grey band show the mean and standard deviation of 1000 simulations 
drawn from an infinite flat model. The WMAP data are consistent with the infinite model. 
The dashed line shows the mean hkelihood, marginahzed over orientation. For the largest 
topology scale, there is no difference between maximizing or marginalizing over orientation; 
the scale is too large for the topology to break the global isotropy. At smaller scales, ori- 
entation starts to matter and although a a parent population may lack global isotropy, any 
given realization may "appear" to break it. The resulting bias can be quantified using Monte 
Carlo simulations. 
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Fig. 4. — Probability Distributions from Monte Carlo simulations, (a) Probability to obtain 
maximum-likelihood output Lgut as a function of the input scale size Lifi after maximizing 
over orientation. The solid line shows Lout = Lin- For Lj„ < 2, the fundamental domain fits 
inside the decoupling surface and the distributions are strongly peaked at Lout — Lin- For 
Lin > 2, maximizing over orientation causes all the simulations to fall into a group centered 
around Lgut ~ 2.2. (b) By inverting the relationship in (a), we obtain the probability that 
an observed output Lout was drawn from a parent population with Lj„. Values Lout > 1-96 
are indistinguishable from infinite models. 
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Fig. 5. — Cumulative probability that the true topology scale is below a given value. WMAP 
data for galaxy cut |6| > 5°. We obtain limits L > 2.1 the distance to the decoupling surface 
at 68% confidence, and L > 1.2 at 95% confidence. The WMAP data is consistent with a 
lack of any finite topology for the cubic fiat topologies we consider. 



